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We present a calculation of the two-loop bottom-sbottom-gluino contributions to Higgs boson 

^\i . production via gluon fusion in the MSSM. The calculation is based on an asymptotic expansion in 

^ I the masses of the supersymmetric particles, which are assumed to be much heavier than the bottom 

\^ I quark and the Higgs bosons. We obtain explicit analytic results that allow for a straightforward 

■^ ■ identification of the dominant contributions in the NLO bottom corrections. We emphasize the in- 

terplay between the calculations of the masses and the production cross sections of the Higgs bosons, 
discussing sensible choices of renormalization scheme for the parameters in the bottom/sbottom 
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1 Introduction 

With the commg into operation of the Large Hadron Colhder (LHC), a new era has begun in the search 
for the Higgs boson(s). At the LHC the main production mechanism for the Standard Model (SM) 
Higgs boson, Hsm, is the loop- induced gluon fusion mechanism [I], 95 — )• -f^sM, where the coupling of the 
gluons to the Higgs is mediated by loops of colored fermions, primarily the top quark. The knowledge 
of this process in the SM includes the full next-to-leading order (NLO) QCD corrections O [3], the 
next-to-next-to- leading order (NNLO) QCD corrections [4] including finite top mass effects [5], soft- 
gluon resummation effects [6|, an estimate of the next-to-next-to- next-to-leading order (NNNLO) QCD 
effects [7| and also the first-order electroweak corrections [51 [9l [10] . 

The Minimal Supersymmetric extension of the Standard Model, or MSSM, features a richer Higgs 
spectrum which consists of two neutral CP-even bosons h, H, one neutral CP-odd boson A and two 
charged scalars H^ . The gluon- fusion process is one of the most important production mechanisms for 
the neutral Higgs bosons, whose couplings to the gluons are mediated by colored fermions and their 
supersymmetric partners. The gluon- fusion cross section in the MSSM is known at the NLO. The 
contributions arising from diagrams with squarks and gluons were first computed under the assumption 
of vanishing Higgs mass in ref. [11]. The complete top/stop contributions, including stop mixing and 
gluino effects, were computed under the same assumption in ref. [12], and the result was cast in a 
compact analytic form in ref. |13) . Later, more refined calculations aimed at the inclusion of the full 
Higgs-mass dependence. In particular, the full squark-gluon contribution is known in a closed analytic 
form |14t ITB| [TB] , while the full quark-squark-gluino contribution has been computed in ref. [T7] via a 
combination of analytic and numerical methods. 

It should be stressed that, at least for the case of the light Higgs, the exact two-loop QCD Higgs- 
gluon-gluon amplitude is in general well approximated by the amplitude evaluated in the limit of 
neglecting the Higgs mass. The latter is much easier to compute and the corresponding result can be 
straightforwardly implemented in computer codes that aim to evaluate the Higgs boson production 
cross section in a fast and efficient way. Indeed, it was noticed several years ago for the SM case [18] 
that the exact K factor, defined as the ratio between the NLO and leading-order (LO) cross sections, 
is well approximated by the so-called effective K factor that can be obtained via an improved effective- 
theory calculation. By the latter we mean a result in which the effective NLO cross section is obtained 
by multiplying the exact LO partonic cross section by the 0{as) corrections evaluated in the limit 
of vanishing Higgs mass. For the SM case this approximation works at the level of few per cent for 
Higgs mass values below the 2mt threshold, and up to 10% for any Higgs mass value. The same level 
of accuracy is reached when the Higgs couples to a generic scalar particle with mass nis, with the 
exception of a narrow region close to the m^ ~ 2 ms threshold [19] . 

There is only one case in which the effective approximation does not work sufficiently well, namely 
when the bottom contribution becomes very relevant. This can happen in the MSSM when tan/3, 
i.e. the ratio of the vacuum expectation values (vev) of the neutral components of the two Higgs 
doublets, becomes large. In such a situation, in principle, the exact computation of the NLO bottom 



contribution to the cross section should be employed. However, explicit analytic results for the part 
of the NLO bottom contribution that is induced by two-loop bottom-sbottom-gluino diagrams have 
not been made available so far. 

In this paper we present an approximate evaluation of the bottom-sbottom-gluino diagrams, based 
on an asymptotic expansion in the large super symmetric masses that is valid up to and including terms 
of 0{ml/mV), 0{mh/M) and 0(m|/M^), where m^ denotes a Higgs-boson mass and M denotes 
a generic superparticle mass (M = rrig^m-j^ ,m-^ ). Together with the known exact results for the 
(s)bottom-gluon diagrams [3l [HI [15l [HI [T9|, our new result allows us to obtain effective K factors 
that can be easily implemented in computer codes to provide an accurate and efficient evaluation of 
the cross section for Higgs boson production in the MSSM. Such K factors are expected to be at the 
same level of accuracy as in the SM, i.e. within a few per cent of the exact value for the light Higgs 
and up to ten per cent for the heavy Higgs, with however a much better accuracy expected if the 
heavy-Higgs mass is below all the heavy-particle thresholds. 

The paper is organized as follows: in section[2]we summarize general results on the cross section for 
Higgs boson production via gluon fusion. In section [3] we outline the way we perform the asymptotic 
expansion of the bottom-sbottom-gluino diagrams. Section H] contains the explicit results for the NLO 
bottom contribution both in DR and in the on-shell (OS) scheme. In section[5]we discuss the numerical 
relevance of the NLO bottom corrections devoting particular attention to the interplay between the 
computation of the Higgs mass and that of its production cross section. In the last section we discuss 
an approximate way to take the NLO bottom contribution into account, based on an improved LO 
term. Finally, in appendix A we specialize to the MSSM case the general exact results for the real 
radiation derived in ref. |19j, while appendix B contains the explicit expressions for the shifts from 
the DR to the OS parameters in the sbottom sector. 

2 Higgs boson production via gluon fusion at NLO in the MSSM 

In this section we recall for completeness some general results on Higgs boson production via gluon 
fusion. The hadronic cross section for Higgs boson production at center-of-mass energy ^/s can be 
written as 

a{hi + h2 ^ (f) + X) = y2 dxidx2fa,hA^i^f^F)fb,h2ix2,fJ-F)x dzSiz ^]aab{z), (1) 

^Jo Jo \ X1X2J 

where (j) = {h,H), r^ = rn^/s, fip is the factorization scale, fa,hi{x,^p) the parton density of the 
colliding hadron hi for the parton of type a (for a = g,q,q), and dab the cross section for the partonic 
subprocess ab ^- ip + X at the center-of-mass energy s = X1X2S = m\/ z. The latter can be written 
in terms of the LO contribution a^^> and a coefficient function Gab{z) as 

aab{z) = a^''^zGab{z). (2) 



We consider now the production of the hghtest CP-even Higgs boson, /i, through gluon fusioru, 
The LO term can be written as 
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where G^ is the muon decay constant, as^^in) is the strong gauge couphng expressed in the MS 
renormahzation scheme at the scale ^u^, Tp = 1/2 is a color factor, and a is the mixing angle in 
the CP-even Higgs sector of the MSSM. T-Li (i = 1,2) are the form factors for the coupling of the 
neutral, CP-even component of the Higgs doublet Hi with two gluons, which we decompose in one- 
and two-loop parts as 
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The one-loop form factors TiY and ?^2^ contain contributions from diagrams involving quarks or 
squarks. The two-loop form factors Til^ and ?^2^ contain contributions from diagrams involving quarks, 
squarks, gluons and gluinos. Focusing on the contributions involving the third-generation quarks and 
squarks, and exploiting the structure of the Higgs-quark-quark and Higgs-squark-squark couplings, 
the form factors Tii can be written to all orders in the strong interactions as [13] 
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In the equations above A( = l/sin/5 and A^ = l/cos/3, where tan/3 = f2/^i is the ratio of the vev of 
the two Higgs doublets. Also, /i is the higgsino mass parameter in the MSSM superpotential, Aq (for 
q = t,b) are the soft SUSY-breaking Higgs-squark-squark couplings and 9q are the left-right squark 
mixing angles (here and thereafter we use the notation s^p = sine/?, c^p = cos<p for a generic angle ip). 
The functions Fq and Gq appearing in eqs. ([5]) and ([6]) denote the contributions controlled by the 
third-generation Yukawa couplings, while Dq denotes the contribution controlled by the electroweak, 
D-term-induced Higgs-squark-squark couplings. The latter can be decomposed as 
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where I^q denotes the third component of the electroweak isospin of the quark q, Qq is the electric 
charge and 9w is the Weinberg angle. 

The one-loop functions entering Til and 7^2 S''^^- 
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^ For the heaviest eigenstate, H, general formulae for the production cross section can be obtained straightforwardly 
with the replacements (sin a — >■ —cos a, cos a — !> sin a) in eqs. (|3]) and (|17|l and in appendix A. 
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where r^ = 4mymf^, and the functions Qq and Gl/2 read 
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The analytic continuations are obtained with the replacement m^ — )■ m| + ie . For later convenience, 
we recall the behavior of Qjf' and ^jL i'^ the limit in which the Higgs boson mass is much smaller or 
much larger than the mass of the particle running in the loop. In the first case, i.e. r >> 1, which 
applies to the top and squark contributions for the light-Higgs case, 
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while in the opposite case, i.e. r <C 1, which is relevant for the bottom quark. 

The coefficient function Gab{z) in eq. ^ can be decomposed, up to NLO terms, as 

Ga,{z) = G^^iz) + ^G^^iz) + 0{al), 

vr 

with the LO contribution given only by the gluon-fusion channel: 
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The NLO terms include, besides the gg channel, also the one-loop induced processes gq 
qq -^ gh: 
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where the LO Altarelli-Parisi splitting functions are 
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In the equations above, Ca = N^ and Cp = (-/Vc ~ l)/(2-^c) [^c being the number of colors), 
/3o = (11 Ca — 2 Nf)/6 {Nf being the number of active flavors) is the one-loop /3-function of the strong 
coupling in the SM, and 

P.W = [i^^iii^l . (20) 

The gg-ch&nnel contribution, eq. (I17p . involves two-loop virtual corrections to gg ^ h and one- 
loop real corrections from gg ^- hg. The former, regularized by the infrared-singular part of the real 
emission cross section, are displayed in the first line of eq. (|17p . The second line contains the non- 
singular contribution from the real gluon emission in the gluon fusion process. The latter contribution 
as well as the ones due to the qq — )• hg annihilation channel and the quark-gluon scattering channel, 
eq. (|18p . are obtained from one-loop diagrams where only quarks or squarks circulate in the loop. 
General expressions for the functions TZgg, Tlqq-, Tlqg can be found in ref. [19] (see also refs. [201 I^T]). 
In appendix A we provide expressions in which the contribution of the bottom quark is kept exact 
while those of the top quark and of the squarks are evaluated in the limit of vanishing Higgs mass. 

The two-loop top/stop contributions to the form factors T-LI^ entering eq. p!7|) are fully under 
control in the light-Higgs case. Typically, the mass ratios between the Higgs and the particles running 
in the loops allow for the evaluation of the relevant diagrams via a Taylor expansion in the Higgs 
mass, with the zero-order term in the series already a very good approximation of the full result. The 
case of the two-loop bottom/sbottom contributions is obviously different. In general, Taylor-expanded 
evaluations of the relevant diagrams are no longer viable, due to the presence of the light bottom quark 
in the loops. Thus, the diagrams must be evaluated either exactly or via an asymptotic expansion in 
a large mass or momentum. 

In the following sections we present the result for the NLO bottom contribution, combining earlier 
results in the literature with our new calculation of the bottom-sbottom-gluino contribution. The 
latter has been obtained via an asymptotic expansion, retaining terms of 0{m^/m^, 0{mi,/M) and 
0(m|/M2). 

3 Outline of the calculation 

An exact analytic evaluation of the bottom-sbottom-gluino contribution to Til 2 is, at the moment, 
beyond our computational ability. However, it is reasonable to assume that all of the supersymmetric 
particles are much heavier than the lightest Higgs boson and the bottom quark, and look for an 
approximate evaluation of the diagrams in terms of a small-momentum (large-mass) expansion. We 
follow this path by performing a large-mass expansion, assuming all the supersymmetric particles to 
be heavy but without requiring any specific hierarchy among them. 

After generating the two-loop diagrams involving bottom, sbottom and gluino that contribute to 
the process g{qi) + 5(^2) — ^ h{q) with the help of FeynArts [22], we separate them in two classes: i) 
those that can be evaluated via an ordinary Taylor expansion in powers of q'^/M'^, of which we keep 
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Figure 1: Pictorial example of the asymptotic expansion of a two- loop diagram containing one subin- 
tegration (ki) with only heavy particles (bold lines) and the other (^2) with light particles. To the 
original diagram its IR-divergent part, represented by the disconnected diagram, is subtracted and 
added forming the contribution that can be evaluated via a Taylor expansion (A) and the one that 
should evaluated exactly (B). See text for a detailed explanation. 



only the term of order zero; ii) the diagrams that require an asymptotic expansion. We recall that 
a Taylor expansion of a two-loop diagram in the external momentum q'^ is viable for values of q'^ up 
to the first physical threshold. In our case, diagrams with a physical threshold at q^ = 47n^, when 
Taylor-expanded in g^, exhibit an infrared (IR) divergent behavior as the bottom mass is sent to zero. 
Thus, these diagrams belong to the class ii. 

Class-i diagrams are expressed in terms of two-loop vacuum integrals that can be evaluated using 
the results of ref. [23]. Concerning the diagrams belonging to class ii, reviews of the method of 
asymptotic expansions of Feynman diagrams with respect to masses or momenta can be found in 
ref. [23]. In practice, we generate the expansion of a diagram by adding and subtracting to it the 
part of the diagram itself that becomes IR-divergent when rrih and q'^ are sent to zero. Formally we 
are adding nothing to the original diagram but, as graphicalljo exemplified in fig. (H this construction 
allows us to separate the diagram in two parts: part A in fig. [Tj which, being by construction IR- 
safe, can be evaluated via a Taylor expansion in the same way as class-i diagrams; part B in fig. [H 
containing the IR-divergent contribution, which should be evaluated exactly. 

The IR-divergent part of a diagram is constructed in the following way. We first note that in all 
the diagrams entering our calculation one can choose a routing of momenta such that the connecting 
propagators, i.e. the propagators that contain both integration momenta fci and /c2, are always accom- 
panied by a heavy mass M. Furthermore, only one subintegration, let us assume the one on k2, is IR 



^The diagrams have been drawn using JaxoDraw [25) 



divergent. Then, one can rewrite the connecting propagators using the identity 

1 _ 1 k1 + 2ki-k2 

(fci + ^2)2 - M2 ~ kf - M2 " [{ki + k2Y - M^]{kf - M2) 



(21) 



The first term on the r.h.s. of eq. (j2ip leads to a disconnected integral (product of two one-loop 
integrals) that contains the IR-divergent contributions present in the original diagram. This term can 
be evaluated exactly, i.e. for arbitrary q'^, giving rise to the ln{q'^ /ml) terms that describe the physical 
threshold. The second term, instead, leads to a two-loop integral with improved infrared convergence 
in the k2 integration and improved ultraviolet convergence in the fci integration. Therefore, if, for 
example, the original integral is logarithmically IR divergent in the k2 integration when g — t- and 
nib —5- 0, the corresponding two-loop integral associated with the second term in eq. (f2T]) evaluated at 
q^ = nib = is no longer IR divergent, but it actually gives a finite result that differs from the result 
valid for q'^ ^ and nib 7^ by terms of 0{m? /M'^ \'a.{m? / M'^)) ^ where m^ denotes either q^ or m^. In 
general, a repeated application of eq. (j2ip . controlled by the power counting in the IR-divergent terms, 
allows us to construct the IR-divergent part of any diagram in terms of products of one-loop integrals 
with numerators that contain terms of the form [ki ■ qj)"^, (ki ■ kj)" {i,j = 1, 2) where ni, n are generic 
powers. The Passarino-Veltman reduction method is then applied to eliminate the numerators and 
express the result in terms of the known one- loop scalar integrals [26] . A check of the validity of our 
construction of the IR-divergent part of a diagram is given by the evaluation of its part A. Indeed, 
one verifies explicitly that the IR-divergent contributions of the original diagram are canceled by the 
terms constructed via eq. (j21|) . so that the final result for part A is free of any h\{q^ /ml^) or q^ /ml 
term. 

4 Two-loop bottom/sbottoni contributions 

In this section we present the result for the two-loop bottom/sbottom contribution to the form factor 
for Higgs boson production via gluon fusion. We stress that, in the MSSM, the result for the production 
cross section of a Higgs boson is strictly linked to the computation of its mass, i.e., both observables 
should be computed in terms of the same set of SUSY parameters, defined in the same way beyond tree 
level. There is however an important difference between the two calculations. In the computation of 
the one-loop corrections to the Higgs masses, the diagrams involving the bottom quark are suppressed 
by the bottom mass and can be safely neglected, resulting in a one-loop contribution that is actually 
due to the sbottom diagrams only. This implies that in the two-loop calculation of the Higgs masses 
the only couplings that require a one-loop renormalization are the trilinear sbottom-Higgs couplings, 
while the definition of the bottom-Higgs Yukawa coupling beyond tree level is irrelevant. On the 
other hand, in the one-loop calculation of the amplitude for Higgs production both the bottom-Higgs 
Yukawa coupling and the trilinear sbottom-Higgs couplings play a role, thus they both require a 
one-loop definition when the two-loop contributions are computed. 
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Figure 2: Examples of two-loop diagrams involving the Higgs-bottom coupling. 
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Figure 3: Examples of two-loop diagrams involving the Higgs-sbottom coupling. 



In the following we will discuss separately the contributions to the two-loop form factors 'H\^2 
from the diagrams involving the Higgs-bottom coupling (examples of which are shown in fig. [2]) and 
those from the diagrams involving the Higgs-sbottom couplings (see fig. [3|). In our approximation 
of neglecting terms beyond 0{mh/M), the former contribute only to the function Gb in eq. ([5]), 
while the latter contribute only to the functions Ff,, Ff, and Gh in eqs. dSJ and ^. For both kinds 
of contributions, we first report the results obtained in the DR renormalization scheme, which is the 
scheme employed by several public computer codes |27j that return the MSSM mass spectrum starting 
from a set of high-energy boundary conditions for the SUSY-breaking parameters. We then discuss 
how to move from the DR scheme to a different renormalization prescription, which generalizes the one 
introduced for the calculation of the Higgs masses in ref. [28] and implemented in the code FeynHiggs 
|29j . With a slight abuse of language, since some of the amplitudes involved in the renormalization 
of the SUSY parameters are in fact evaluated off mass shell, we refer to this prescription as to the 
"on-shell" scheme. 



4.1 Contributions controlled by the Higgs-bottom coupling 

We start by considering the contributions of the two-loop diagrams involving the Higgs-bottom cou- 
pling. In our approximation those are the only diagrams that contribute to the function Gb in eq. ([5]). 
The two-loop part of the function can be decomposed as 
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Assuming that the one-loop form factor T-LY is expressed in terms of DR-renormalized parameters 
evaluated at the scale Q^, the contribution of the two- loop diagrams with bottom quarks and gluons 



(fig. H^) reads: 



2 r^{g,CF) 



2miG, 



-^ir\-5) + -^ir(-.) In 



n 



ml 1 
Q2 "3 



2miG 



2 r<i9,CA) 



{2i,CA)i 



lGf'--> = gY;^-\n). (23) 



Exact expressions for the functions ^112 ■, -^1/2' ^^^ ^1/2 ^ ^'^^ given in eqs. (2.12), (2.13) and 
(3.8) of ref. [K], respectively. In the hmit r ^ 1 they reduce to 
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where C2 and Cs are Riemann's zeta functions. 

The contributions of the two-loop diagrams with bottom, sbottom and gluino (fig. [2}3) require a 
dedicated calculation. Up to and including terms of 0(m^/?n,^) and 0{mh/M), and assuming that 
TiY is expressed in terms of DR-renormalized parameters evaluated at the scale Q^, they read: 
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where Xj = m? /m| , and {6mb)^^^^ denotes the SUSY contribution to the bottom self-energy, in units 
of Cp Os/tt and in the limit of vanishing rrib- 
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where 



If the bottom-quark contribution to Tii is expressed in terms of the pole bottom mass Mb, the 
functions G^^' ^ and G^^' ^ are shifted with respect to their expressions in eqs. (|23p and (|27|) . In 
particular, the former becomes 



2mlGr'^^ = ^(J")(,,) + ^^(J^)(,,), (31) 



and the term proportional to {dTUb)^^^'^ in the first line of eq. (f27]) is canceled out. 

Eqs. (|27p and (j29p show that the bottom-sbottom-gluino contribution to Til contains terms en- 
hanced by the large ratio rrig/mf,. Recalling the definition of r^, it is clear that those terms are in fact 
of 0{mbmg/m\) , i.e., they still vanish as mj, — )• but they are enhanced by the ratio mg/rrih. Such 
terms arise from two-loop diagrams in which the helicity flip on the fermion loop is achieved via a 
gluino mass insertion instead of a bottom mass insertion, and they by far dominate the new-physics 
contribution to the two-loop part of the form factors, no matter whether the bottom-quark contribu- 
tion to I-LY is expressed in terms of the pole bottom mass M^ or in terms of the DR-renormalized 
bottom mass "fhb- However, we notice that all of the two-loop 0{mbmg/mf^) terms cancel out if the 
one- loop bottom contribution to the function G^^ is computed in terms of Mb, but the function itself 
is multiplied by ffibMb instead of M^. As a result the function G^^ is further shifted, with respect to 
the expression corresponding to the use of Mb in the whole one-loop contribution, by 
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This manipulation amounts to differentiating, in the one-loop contribution, between the parameter 
that describes the mass of the bottom quark running in the loop - which is identified with Mb - and 
the parameter that describes the Yukawa coupling of the bottom quark to the Higgs boson - which is 
identified with mj,. We recall that, in the MSSM, the running bottom mass irib can be related to the 
corresponding SM parameter m^ as [30] 

mb{l + db) , . 

^b = :r- ; o > (33) 

I + eb tan p 

where 5b denotes terms that are not enhanced by tan/3 and, to 0{as), 

asGp 2fimg / xi X2 \ , . 

eb = —. 2 2 ^i li^ ^1 ~ ~i In X2 . (34) 

47r mi — m^ \1 — xi 1 — X2 / 

01 02 

Since S2e^ = 2mb {Ab + /itan/3)/(?7i? — m? ) , it is easy to see that the terms enhanced by m.g/m,b in 

Gff' ^ do indeed contain eb tan/3. In the effective-theory language of ref. [31] we can argue that, by 
expressing the bottom Yukawa coupling entering TiY in terms of fhb as defined in eq. ()33p . we "resum" 
in the one- loop part of the form factor the tan /3-enhanced threshold corrections to the relation between 
the mass and the Yukawa coupling of the bottom quark. As a result of this special choice of parameters, 
all terms of the form eb tan /3 drop out of the two-loop part of the form factor. 
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4.2 Contributions controlled by the Higgs-sbottom coupling 

We now turn our attention to the diagrams that involve the Higgs-sbottom couphng. In our approxi- 
mation those diagrams contribute only to the functions i*),, i*), and G^ in eqs. ([5]) and ([6|). In analogy 
with ref. [13j, the two-loop parts of the functions can be written as 

^f = n, - Yi. + J'"'' 2 ya. > (36) 



61 62 



Gf = n,+n2- (37) 

Furthermore, the various terms in eqs. ()35p - (|37p can be split in the contributions coming from diagrams 
with gluons [g, fig. [3^), with strong, D-term-induced quartic bottom couplings (46, fig. [3)3), and with 
gluinos {g, fig. Et), 

y, = Y3 + Y^^ + Yi {x = h,b2, cL) . (38) 



The first two terms in eq. (|38p can be obtained from the analytic expressions presented for the stop 
contributions in eqs. (27)-(30) of ref. [13], identifying ^- with Y^ after making the trivial replacement 
t ^ b. The gluino contributions, on the other hand, require a dedicated calculation. Writing 

Yi = Cp Yf3,c,) ^ Ca yj^'^^) {x = 61, 62, c2,J , (39) 

and assuming that the parameters in T-Li and T-O^ are expressed in the DR scheme at the renormal- 
ization scale Q^, we find for the functions Y^: 

^t"'^ = r^^^i/2(-^)fr^ + 7r^i--i 

''1 Ambuig '-'^ Vl-2;i (1 - xi)^ 

1/1 1 11, mi\ , , 

+ JZ -^Inxi-^ + ^ln-f , (40) 



6m? ll — xi (1 — xi)2 xf xf Q"^ 



9 
1/1 1 



^1 12m?\l — xi (1 — xi)2 



n'*''''' = -TT^Ir^^ + TT^rrjta^il . («) 



fF" = -^^=L5;;^(,,)(_S_i,„,-^,„.,), (42) 

"^Ob 8'mbS2et ' Vl-a:;i 1 - X2 / 

y|'^^) = 0, (43) 

where we retained only terms that induce ©(m^/ml), 0{mb/M) and 0{m'^/M'^) contributions to 
T-Li2- The expression for Yj^^ can be obtained from the expression for Yj^ through the replacements 
xi — )• X2 and 829^ — ^ "■5206- Comparing eqs. (|iO|) and ([^T]) with eq. (42) of ref. |13] we notice that, 
contrary to what we stated in section 3.3 of that paper, even for Of, = the two-loop bottom-sbottom- 
gluino contribution to Fj, and Gj, cannot be obtained by taking the limit mj — )• in the corresponding 
top contribution. 
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4.3 On-shell renormalization scheme for the sbottom parameters 

We now discuss a suitable OS renormalization scheme for the parameters that determine the sbottom 
contribution to ^{^ We recall that, at the one-loop level, the vev vi and V2, the Z boson mass, the 
Weinberg angle and the parameter fj, are not renormalized by the strong interactions. Therefore, the 
only parameters that require a one-loop definition are {hb, Ab, S26»(,j w-^ , m^; ), where by hb we denote 
the coupling constant entering the cubic and quartic sbottom-Higgs interactions, which at tree level is 
related to the bottom mass by rrib = hbVi/ \f2. Indeed, the factor rrib that multiplies the function Fb 
in eqs. ([5|) and ([6]) has to be interpreted as a bookmark for hb- In fact, only four of those parameters 
are independent, because of the relation 

^hb^AbVi^ IIV2) .... 

826^ = 2 2 • (44) 

mf — mf 

In the analysis of the sbottom corrections to the neutral MSSM Higgs boson masses presented in 
ref. [28] it was pointed out that, while the sbottom masses can be naturally identified with the pole 
masses, an OS definition of {hb,Ab,S2ei,) is less easily singled out. Proceeding in analogy with the OS 
renormalization of the stop sector (see, e.g., ref. |13]). we might choose as independent parameters 
a conveniently defined bottom mixing angle, 5261;,) f^d the bottom Yukawa coupling /i^"*^, as defined 
by the pole bottom mass Mb via the relation Mb = h^°'^vi/V^. Then, eq. (fill) might be used to 
establish the one-loop definition of Ab in terms of the pole bottom and sbottom masses and the 
sbottom mixing angle. However, for large values of tan/3 such definition would produce very large 
shifts in Ab with respect to its DR valueQ 6Ab = 0{as^^ i&v? j3/mg) [32]. This is related to the 
fact that, in the large-tan/3 limit (i.e., vi — )• 0), 829,, becomes independent of Ab^ as can be easily 
seen from eq. (HH). To cure the problem, it was suggested in ref. [28] (see also ref. |33J) to take 526^ 
and Ab as independent parameters, while considering hb as a derived quantity via eq. (jM]). Suitable 
renormalization conditions were then proposed for S261,, and Ab- 

In the OS analysis of the cross section for Higgs boson production we want to retain the convenient 
features of the renormalization prescription employed in ref. [28j. However, that prescription needs to 
be expanded: first of all, the renormalization conditions in ref. p8] were defined in the limit tan /3 — )■ 00 
(i.e., vi =0), while in the case at hand we do not impose constraints on tan/3. Moreover, while in the 
calculation of the one-loop corrections to the Higgs masses the contributions controlled by the bottom- 
Higgs Yukawa coupling (which we denote as }i( to distinguish it from hb) are suppressed, the one-loop 
diagram controlled by h^ gives an important contribution to the production cross section. Therefore, 
a one-loop definition of h^ is required. Since it does not seem appropriate to define the bottom-Higgs 
coupling /i^ in terms of quantities of the sbottom sector, as would happen if we imposed on it the 
same renormalization condition used for hb, we need to impose different renormalization conditions on 
hi and hb, or, equivalently, on the bottom mass that enters the one- loop bottom contribution and the 
one that enters the one-loop sbottom contribution. In particular, we identify the former with the pole 



^For the generic parameter x, we define tiie shift from the DR value x as 5x ^ x — x. 
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mass Ml, (the resulting shifts in the function G^ are discussed in section H7T]) . Concerning the bottom 
mass in the one-loop sbottom contribution, we follow ref. [28], extending the prescription presented 
in that paper to the case of finite tan/3. 

To obtain definitions for J/ij and SA^, we consider two quantities 

Xb = —i=- {Ab + iJ, tan P) , % = — (si3 Ab - cp n) , (45) 

that allow for a natural interpretation: Xb , at the classical level, is the off-diagonal term in the sbottom 
mass matrix, related to the mixing angle S2e^ via eq. (|44p ; Yj, is proportional to the coefficient of the 
trilinear interaction (6162^). A definition of the mixing angle 6b like the one proposed in ref. [34.] , 

1 ni2(m?) + ni2(m? ) 

Seb = i ^ ^^^ , (46) 

2 mf — mf 

together with the identification of the sbottom masses as pole masses, can be immediately translated, 
using eq. ([Bl), into a prescription for Xb'- 



SXb ^ 



1 r^ . o. ^ . 0.1 -~, Uu{ml)-U22{ml 



b2' 



b2' 



ni2(m )+ni2(m J +Xb ^i -2 ^. (47) 



2 2 



bl 62 

In eqs. (jl6|) and ([T7|) . IIjj {q'^) (i, j = 1, 2) denotes the finite part of the {i,j) self-energy of the sbottoms. 
Recalling that in the tan/3 — )• 00 limit Yb — ^ hbAb/V2, the extension to the case of finite tan/3 of 
the prescription for Ab introduced in eq. (15) of ref. [2Sj reads: 

SYb = - - [Ai2A (ml^ ,ml^,0)+Ai2A (m?^ , m?^ , 0) 

1 ^ nii(m? ) + U22{ml ) - nii(m? ) - 1122 (m? ) 

+-Yb '-^ % -, '-2 ^ (48) 

2 mf — m,f 

where iAi2A{Pi, P2j Pa) denotes the proper vertex 6162^. 

Finally, the shifts of the parameters hb and Ab are related to those of Xb and Yj, by 

5hb = —isXbSp-SYbVCp) , (49) 

^^b = -^(Xb6Yb-Yb6Xb) , (50) 



where v = Jv\ + t;|. Explicit expressions for 5hb and 5Ab, as well as for 5s2ef^ and (5rM? , can be found 
in appendix B. 

If the one-loop sbottom contribution to l-L^ ^^^ ^2^ is evaluated in terms of OS quantities, the 
two-loop functions in eqs. (j35]) - (j37j) must be replaced by 
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(51) 



(52) 



(53) 



In addition, the two-loop form factor T-Lf receives a contribution originating from the shift in A},: 



nf 



nf 






vr 
6 a., 



rrir 



rrir 



SAb 



(54) 
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5 A numerical example 



We win now illustrate the effect of the two-loop bottom/sbottom contributions to the form factors for 
the production of a Higgs boson in a representative region of the MSSM parameter space. 

The SM parameters entering our calculation are the Z boson mass niz = 91.1876 GeV, the Fermi 
parameter Gp = 1.16637 x 10~^ GeV~^, the sine of the Weinberg angle Sq = 0.223 and the strong 
coupling constant Osimz) = 0.118 [35]. For the pole masses of the top and bottom quarks we take 
Mt = 173.1 GeV [36j and Mb = 4.49 GeV, the latter corresponding to the SM running mass (in the 
MS scheme) rnb{mb) = 4.16 GeV |37] . The tree-level mass matrix for the CP-even Higgs bosons can be 
expressed in terms of the physical pseudoscalar mass tjia and the DR-renormalized parameter tan/3, 
in addition to mz- In the calculation of the physical Higgs boson masses and of the mixing angle a 
we include the one-loop 0{at + a^) and two-loop 0{atas + abOig) corrections as in refs. [251 [35] • 

When computing the two-loop corrections to both mass matrix and production form factors for the 
CP-even Higgs bosons, the parameters that determine the stop and sbottom masses and mixing angle 
and are subject to 0{as) corrections require a one-loop specification. For the stop sector we adopt 
the OS scheme described e.g. in ref. [T3]. In particular, we take as input the pole top mass Mt and the 
soft SUSY-breaking parameters {mQi,mu,At) that can be derived by rotating the diagonal matrix 
of the OS stop masses by the angle Ot, defined as in eq. (37) of ref. [13]. Concerning the corresponding 
parameters of sbottom sector {hb, tUq ^ , mj:i,Ab) additional care is required, because of our non-trivial 
definition of hb and of the fact that, at 0{as), the parameter rriQ ^ entering the sbottom mass matrix 
differs from the corresponding stop parameter ran t by a finite shift |32j. We start by computing the 
renormalized Higgs-sbottom coupling as given by hb = hb — 5hb, where hb is the DR-renormalized 
running coupling that can be trivially extracted from irib computed via eq. ([33]) . and 5hb is defined in 
eq. (j49|) . Then we compute tjIq ^ following the prescription of [32] . Finally, we use the parameters hb 
and mp ^ to compute the actual values of the OS sbottom masses and mixing angle. 
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Figure 4: CP-even Higgs boson masses (left plot) and effective couplings of h to top and bottom 
quarks (right plot) as a function of tan/3, for mA = 150 GeV and all SUSY mass parameters equal to 
M = 500 GeV. For the meaning of the solid and dotted lines see the text. 



To set the stage for further discussion, we show in fig. [J] the masses and mixing angle of the 
CP-even Higgs bosons as a function of tan/3. All the relevant SUSY-breaking parameters, as well 
as the supersymmetric mass parameter /i, are set to a common value M = 500 GeV, and the phys- 
ical pseudoscalar mass ruA is set to 150 GeV. The left panel of fig. H] shows the masses of the two 
Higgs bosons h and H in the range 25 < tan/3 < 50, while the right panel shows the combinations 
— sin a/ cos /3 and cosa/sin/3, which determine the strength of the coupling of h to the bottom and 
top quarks, respectively, relative to the corresponding SM couplings. For each set of curves, the solid 
line represents the result obtained in the OS renormalization scheme of ref. [28], described in section 
14. 3[ For comparison, we also show as a dotted line the result that would be obtained if the sbottom 
parameters hb and Af, were renormalized in the same way as the corresponding stop parameters. The 
left plot shows the well-known fact that, at large tan/3, the radiative corrections from sbottom loops 
tend to reduce m^. The right plot shows that, for the chosen values of rriA and tan/3, the coupling 
of h to the bottom quark is still substantially enhanced with respect to its SM value. This has to be 
contrasted with the couplings of h to the top quark and to the gauge bosons (not shown), which are 
already very close to the SM values they tend to in the "decoupling" limit ruA ^ rriz- 

The comparison between the dotted and solid lines in fig. H] shows that, if we had adopted for the 
sbottom parameters h^ and Af, the renormalization scheme used for the stop parameters, the results for 
niH and for — sin a/ cos /3 would differ wildly from the ones obtained with the renormalization scheme 
discussed in section [3]3] (conversely, we checked that the results obtained in the DR scheme would be in 
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Figure 5: K factor for the production of a light Higgs boson h as a function of tan/3, for rriA = 150 
GeV and ah SUSY mass parameters equal to M = 500 GeV. For the meaning of the different lines 
see the text. 



good qualitative agreement with the solid lines). The discrepancy is due to the fact that in the "bad" 
OS scheme the (1,1) and (1,2) entries of the CP-even Higgs mass matrix are subject to very large 
two-loop corrections scaling like M^tan/3^, induced by the contribution of the counterterm 6Ah. It is 
interesting to note that, since the contribution of 6Ai, to the form factor Til^ in eq. ([5^ is suppressed 
by a factor rub, its impact on the Higgs boson production cross section in the "bad" OS scheme is not 
as extreme as the impact on the Higgs mass. However, we stress that a consistent determination of 
the properties of the Higgs bosons requires that the same definition of input parameters be used in 
the calculations of mass and production cross section. Since the naive choice of using the same OS 
renormalization scheme for the stop and sbottom sectors is not viable in the calculation of the Higgs 
masses, it should not be applied to the calculation of the cross section either. 

We are now ready to discuss the effect of the two-loop bottom/sbottom contributions to the form 
factor for Higgs-boson production. To this purpose, we define a factor Kh that contains the ratio of 
two-loop to one-loop form factors appearing in eq. ([1] 



Kh 



1 + 2 — Re 

n 



sin a T-LY + cos a Ti^^ 



sin aTii + cos a T-L^ 



u 



(55) 



In the left panel of fig. [5] we plot K^ as a function of tan/3, with the same choice of SUSY 
parameters as in fig. [U in the OS renormalization scheme described in section [¥31 The one-loop form 
factors in the denominator of the term between parentheses in eq. (I55p contain both the top/stop and 
bottom/sbottom contributions, computed under the approximations of eqs. (113^ and (I14p . The lines 
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in the plot correspond to different computations of the two-loop form factors in the numerator: the 
dotted line includes only the contributions of the top/stop sector, as computed in ref. [13j; the dashed 
line includes also the contribution of two-loop diagrams with bottom quarks and gluons; the solid 
line includes the full two-loop contribution of the bottom/sbottom sector as computed in section \^ 
finally, the dot-dashed line is obtained by approximating the bottom/sbottom contribution (with the 
exception of the bottom-gluon diagrams) with just the terms enhanced by rrig/rrn, in eq. (I27p . Prom 
the comparison between the dotted and dashed lines it can be seen that, in the OS renormalization 
scheme, the contribution to T-Li of the two-loop diagrams with bottom quarks and gluons is very 
small. This is due to a partial cancellation between the terms CpT^,^ and Ca ^1/2' ^ entering the 
function G^ , and to the fact that, in this scheme, the term T^i^ is not enhanced by the potentially 
large logarithm of the ratio between the bottom mass and the renormalization scale, as can be seen by 
comparing eqs. (j23p and (j3ip . The solid line shows that the effect of the diagrams involving sbottoms 
can be very sizable at large tan/3, more than doubling K^. Indeed, for large tan/3 the coupling of the 
light Higgs boson to the (s)bottom is considerably enhanced with respect to its SM value, as can be 
seen in the right panel of fig. |H However, the proximity between the solid and dot-dashed lines shows 
that this sizable effect is almost entirely due to the terms enhanced by rrig/mh in the contribution of 
the two-loop bottom-sbottom-gluino diagrams in which the light Higgs boson couples to the bottom 
quark. 

As discussed in section 14.11 the terms enhanced by rrig/mi, in the OS result can be canceled out 
if the Higgs-bottom Yukawa coupling in the one-loop part of the result is identified with the DR- 
renormalized MSSM bottom mass rhh instead of the physical mass M^. To this effect, the factor m^ 
multiplying the function G^ in eq. (^ must be expressed as fhbMb, and the two-loop part of Gi, must 
be shifted as in eq. (j32p . In the right panel of fig. [5] we present the result of this manipulation, with fhf, 
evaluated at the scale Q = ruh- The input parameters and the meaning of the different lines are the 
same as for the plot in the left panel. The proximity between the dashed and solid lines shows that 
the contribution of the two-loop diagrams involving sbottoms is rather small in this renormalization 
scheme, at least for our choice of input parameters. However, K^ still shows a sizable increase at large 
tan/3. This is due to the fact that the shift in eq. (|32p brings back a large logarithm, ln(77i^/m|), in 
the contribution of the two-loop diagrams with bottom and gluon (this logarithm compensates the 
scale dependence of the running mass irn,). 

To conclude this section, we show in fig. [6] the factor Kh for the production of the heavy CP-even 
Higgs boson, in the range 2 < tan /3 < 50. The definition of Kh can be obtained from the one of K^ 
in eq. ()55p via the replacements sin a — t- —cos a, cos a — t- sina. The input parameters are chosen 
exactly as in figs. H] and O and the meaning of the lines in the left and right panels is the same as 
in fig. [5j Since in this example the mass of the heavy Higgs boson is of the order of 150 GeV (see 
fig. HI), i.e. well below any threshold for heavy-particle production, we expect the approximation of 
vanishing Higgs mass to hold reasonably well even for H. Prom fig. [6] it appears that the balance of 
the various contributions to Kh in the two different renormalization schemes is qualitatively similar 
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Figure 6: Same as figure [5] for the heavy Higgs boson H. 



to the one for K^ shown in fig. [D in the OS scheme the factor Ku receives a sizable contribution from 
the sbottom diagrams, largely dominated by the terms enhanced by rrig/mi, in the diagrams controlled 
by the Higgs-bottom coupling; in the "mixed" scheme, on the other hand, the sbottom contribution 
is rather small, but there is a sizable contribution from the diagrams with bottom and gluon. 

This said, the factor Kh shows a peculiar dependence on tan/3: for sufficiently large values of 
tan/3, it grows linearly in the OS scheme, while it reaches a plateau in the mixed scheme. This can be 
easily understood by recalling that, for moderate-to-large tan/3 and for our choice of m^, the Yukawa 
coupling of the heavy MSSM Higgs to bottom quarks is enhanced by tan j3 with respect to the SM 
value, while the coupling to top quarks is suppressed by tan /3. Consequently, both the one-loop and 
the two-loop form factors in Kh are dominated by the contribution of the diagrams controlled by 
the Higgs-bottom coupling, with the result that the coupling itself cancels out in the ratio. However, 
the dominant contribution from the bottom-sbottom-gluino diagrams in the OS scheme contains an 
additional tan /3-enhancement hidden in the product S20^'mg/mh (see the discussion at the end of 
section 14. ip . which explains the linear rise of Kh- On the other hand, the dominant contribution 
of the bottom-gluon diagrams in the mixed scheme possesses no further tan /3-enhancement, which 
explains the plateau. 

6 Conclusions and discussion 

In this paper we presented analytic results for the NLO bottom-sbottom-gluino contribution to the 
cross section for Higgs boson production in gluon fusion, obtained using an asymptotic expansion in 
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the large supersymmetric masses. This approximation is fully valid for the light-Higgs case, while 
for the heavy Higgs it covers the mass region where rrin is below all the heavy-particle thresholds. 
Together with the previously known results for the NLO corrections in the MSSM, our expressions 
can be easily implemented in computer codes that aim to provide an accurate and efficient evaluation 
of the cross section for Higgs boson production in the MSSM. 

In our analysis we paid special attention to the consistency between the calculations of the masses 
and the production cross sections of the MSSM Higgs bosons, i.e. to the fact that the same input 
parameters, defined at the one-loop level, should be used in both calculations. The OS definition of 
the parameters of the bottom sector is delicate, as discussed for the case of the Higgs masses in ref. |28j . 
The choice of treating the top and the bottom sectors on the same footing suffers from the fact that 
large two-loop corrections proportional to tan^ /? are generated in the contributions controlled by the 
Higgs-sbottom couplings, affecting both the calculation of the Higgs masses and that of the production 
cross sections. To avoid such large two-loop effects, a convenient OS renormalization prescription was 
proposed in ref. [28] for the calculation of the Higgs masses. In the present paper we have extended 
that prescription to cover also the calculation of the production cross sections. 

Our analysis of the NLO bottom contribution to the gluon-fusion production cross section shows 
that, with our choice of OS renormalization conditions, the bulk of the corrections comes from the 
two-loop diagrams involving the Higgs-bottom Yukawa coupling, while the diagrams controlled by 
the Higgs-sbottom coupling play a secondary role. The contribution controlled by the Higgs-bottom 
Yukawa coupling can be further divided in two parts: diagrams with only bottom and gluons and 
diagrams involving bottom, sbottom and gluino. By far, the most important pieces of the latter 
diagrams are the terms of 0{rrH)mg/mf^ jj), i.e. the ones in which the helicity flip on the fermion line 
is achieved via a gluino mass insertion instead of a bottom mass insertion. 

It is natural to wonder if it is possible to absorb most of the NLO bottom contribution into the LO 
term with a suitable choice of the input parameters. In such a situation the factor K^ij, as defined in 
eq. ([55]) , would be basically sensitive to the top/stop contribution only. The contribution of the two- 
loop bottom-gluon diagrams can be made small if the one-loop bottom diagrams are expressed in terms 
of the pole bottom mass M^, but in this case the 0{mi,mg/mf^ jj) terms give a sizable contribution. 
On the other hand, if the Higgs-bottom Yukawa coupling in the one-loop result is expressed in terms of 
the running bottom mass rhh, and the bottom mass in terms of Mh, the 0{mijmg/m\ ^) terms in the 
two-loop contribution cancel out, but the bottom-gluon diagrams give a relevant contribution because 
of the presence of large logarithms of the ratio between rri}, and the renormalization scale. However, 
the explicit knowledge of the NLO bottom contribution allows us to devise a simple recipe to absorb 
the bulk of the NLO contribution into the LO term. It amounts to writing the LO bottom contribution 
entirely in terms of the pole bottom mass M;,, then rescaling it by a factor 1/(1 -|- e(,tan/3). Once 
this manipulation is implemented, we expect the remaining NLO bottom/sbottom contributions to be 
quite small - at least in large regions of the parameter space - in which case they can be neglected in 
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the evaluation of the form factors ^12 without introducing large errorqj. We stress that the validity of 
this simple recipe is strictly linked to the absence of spuriously large corrections to the Higgs-sbottom 
coupling. This is realized with our choice of OS renormalization conditions for the sbottom sector 
(and also in the DR scheme) but it is not guaranteed with other renormalization conditions. 

Finally, the results derived in this paper for the production cross section can be straightforwardly 
applied to the NLO computation of the gluonic and photonic decay widths of the CP-even Higgs boson 
in the MSSM, as described in section 5 of ref. [13] . 
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Appendix A: NLO contributions from real parton emission 

In this appendix we specialize to the MSSM case the general exact results of ref. [19j for the functions 
TZgg, TZqq, Tlqg- We aim at expressions that, on one hand, are sufficiently accurate, while on the other 
hand allow for a fast numerical evaluation. Thus we report expressions in which the contributions of 
the top quark and of the squarks are evaluated in the limit of neglecting the Higgs mass, while the 
contribution of the bottom quark is kept exact. 



The function TZgg can be written as 



n 



1 



99 



dv 



'^ggx^j '•; ^j 



z{l-z)Jo v{l-v) ] l-sinanY + COS an^f 
-s(l — z)(l — v), u = — s(l — z)v, with 



(1-z- 



.2x2 



(Al) 



where t 

\Agg{s,t,u)\'^ = \A2{s,t,u)\^ + \A2{u,s,t)\^ + \A2{t,u,s)\^ + |A4(s,t,n)|2. 
Furthermore, the functions A2 and ^4 can be cast in the following form: 



A2is,t,u) 
Ai{s,t,u) 



1A2 



Ml 



sin a IZ^ ^{s,t,u) + cos a 7^2 ^ (-^j *> u) , 



sin aTZi"^ {s, t, u) + cos aTl2'^{s,t,u) , 
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*A somewhat similar procedure was suggested, without a detailed discussion, in ref. [39| 
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where the functions 6i/2(-s,i,ii) and Ci/2(s,i,u) are defined in eqs. (2.22) and (2.24) of ref. |T9], respec- 
tively, and it is understood that the top and squark contributions to 'HY2 ^^^ evaluated in the limit 
of neglecting the Higgs mass. In several cases the terms proportional to A;, in eqs. (|A5P and ()A7p are 
numerically very small and can be neglected. In such a situation the integration in eq. ()Aip can be 
performed analytically, resulting in TZgg = —11(1 — z)^/(6z). 
The qq — )• Hg annihilation channel can be written as 

2 



7^, 



128 



zil 



•^qqx^i ^5 ^/ 



qq 



27 |— sina'HJ^ + cosa'H2^| 



with 



Aqq{s, t,u) = — sin Q 7^]^ '' (s, t, u) + cos Q 7^2 '"^ ('^j ^) "") 



where 
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t + u 
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'^br Is 

jdl/2 



2 ' 2 ' 2 

771^ m^ m^ 



+ 



t + u 



2{s + t + u) 



Gl%{n) 



n\' 



2{s + t + u) 

The function dii2{s.,t,u) is defined in eq. (2.31) of ref . [19] . 

Finally, we consider the quark-gluon scattering channel, qg 
be written as 



(A9) 

(AlO) 

, (All) 
(A12) 
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2„,2 



l + {l- zYv 
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•A.qg\S^ t, U^ 
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where 



[l-(l-z)f]2 |_sina-Hj^ + cosa7^^ 



•^qg{S, t, U) — ^qq(t, S, U) 



qH. The relevant function TZqg can 

(A13) 



2z 



(A14) 



As in the case ofTZgg, when the terms proportional to A;, in eq. (|Alip can be neglected the integration 
in eq. (jAlSp can be performed analytically, giving TZqg = 2z/3 — {l — z)'^/z and TZqg = 32(1 — z)'^/(272:). 
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Appendix B: Renormalization scheme shifts in the sbottom sector 



In this appendix we present explicit expressions for the shifts from the DR to the OS scheme of the 
parameters in the sbottom sector that require a one-loop definition. Denoting, generically, a quantity 
in the DR scheme as x°^, and the same quantity in the OS scheme as x°^ , we can write the one-loop 
relation as x°^ = x'^^ + 6x. Retaining only terms that do not induce contributions suppressed by 
ml/M'^, we find: 
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where the notation (1 o 2) in eqs. (jBSp and (|B4p means a term that is obtained from the previous 



ones inside the square bracket with the exchange m? -H- mif . The shift 6m,f is obtained from eq. ( IBll) 

oi 02 02 _^ ^ ^ 

via the interchange m? -f-T- m? . Finally we note that the expressions for 5hb and 6Ab in eq. (jBSp and 



eq. ()B4p . respectively, which are valid for generic values of tan/3, coincide with the corresponding 
expressions in ref. [28], which were derived in the limit tan/3 — )• cxd. 
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